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摘 要 : 


展 矿 区 生态 系统 服务 价值 的 时 空 分 析 对 保护 矿区 生态 环境 、 提 升 矿区 生态 系统 功能 具有 


重要 意义 。 基 于 还 感 数 据 和 土地 利用 数据 ,采用 3S 技 术 、 价 值 当量 法 和 生态 贡献 度 等 方法 ,对 黄土 
高 原 典 型 矿区 一 一 平 朔 矿区 1990 一 2020 年 生态 系统 服务 价值 进行 研究 。 结 果 表 明 :(1) 1990 一 


2020 年 平 朔 矿区 6 种 土地 利用 类 型 面积 变化 角 蜡 ,耕地 


名 减少 最 大 (1122.72 hm?) ,建设 用 地 面积 


增加 最 大 (2044.23 hm?) 。(2) 平 朔 矿区 生态 系统 服务 价值 呈 先 下 降 后 上 升 态 势 , 由 1990 年 的 
19562.43x10+ 元 减少 到 2010 年 的 11265.40x104 元 ,再 增加 到 2020 年 的 13755.47x10: 元 , 1990—2010 
年 下 降 的 主要 原因 是 建设 用 地 和 和 草地 的 生态 系统 服务 价值 大 幅 减 少 ,2010 一 2020 年 增加 的 主要 原 


因 是 矿区 实施 了 大 量 生态 恢复 治理 工程 ,(3) 1990 一 2020 年 平 朔 矿区 除 保持 土壤 价值 增加 外 ,其 他 


各 单项 生态 服务 价值 均 有 不 同 程度 的 减少 ,其 中 水 源 涵养 .净化 环境 、 生 物 多 样 性 价值 减少 幅 
度 较 大 ,分 别 减少 了 115.84%、69.92% 、18.29% 。(4) 1990 一 2020 年 平 朔 矿区 生态 系统 服务 价值 


的 空间 分 布 呈 现 南 北 高 、 中 部 低 的 特征 。(5) 3 


FF 朔 矿 区 草地 、 建 设 用 地 和 林地 的 生态 贡献 度 分别 


是 -41.52% 、-34.49% 和 -10.09% ,三 者 是 主要 的 贡献 因子 和 敏感 因子 。 
X 键 词 : 生态 系统 服务 价值 ， 时 空 变化 ; 平 朔 矿区 ; 黄土 高 原 
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生态 系统 是 人 类 生存 和 繁衍 的 基础 ,为 人 类 提 
供 了 赖 以 发 展 的 环境 条 件 和 物质 基础 "。 生 态 系 统 
服务 是 指 人 类 直接 或 间接 从 生态 系统 的 结构 ,过程 
和 功能 中 获取 的 各 种 产品 和 服务 ,一 般 分 为 供给 、 
支持 .调节 文化 4 类 服务 ”。 生 态 系统 服务 是 联系 
人 类 活动 和 生态 系统 的 桥梁 和 纽带 ,是 研究 人 类 系 
统 -生态 系统 耦合 的 重点 内 容 ”。1970 年 Wilson 首 
次 提出 “生态 系统 服务 ”概念 ,1997 年 Daily ^f 
Costanza 等 5 对 该 概念 进行 了 系统 研究 ,之 后 迅速 引 
起 国内 外 生态 学 、 地 理学 和 环境 学 等 领域 的 专家 学 
者 广泛 关注 ,不 同 空间 尺度 (全 球 ” 国家 "省 市 、 
县 区 外 谷地 高原" 等 ) 和 不 同 生态 系统 类 型 ( 城 
市 中 农田 外 河流 森林 "草原 中 湿地 "等 ) 


收 稿 日 期 : 2023-07-15; ”修订 日 期 : 2023-08-28 


下 的 生态 系统 服务 相关 工作 得 到 普遍 认识 和 研究 ， 
并 取得 较为 丰硕 的 学 术 成 果 , 这 些 成 果 主 要 集中 在 
生态 系统 服务 的 价值 评估 00 ud DA, 
时 空格 局 ”权衡 协同 关系 下、 形成 机 制 "™、 人 
类 福 社 等 方面 。 近 年 来 ,生态 系统 服务 的 研究 由 
自然 生态 系统 (如 :湿地 ) 延 伸 到 半 自 然 生 态 系统 
n: REA), 矿区 是 以 资源 开采 为 核心 的 一 个 
独特 的 半 自 然 半 人 工 生态 系统 ,长 期 资源 开采 造 
成 了 土地 损毁 ,水土 流失 ,水体 污染 植被 退化 等 生 
态 环境 问题 , 极 大 地 影响 了 矿区 生态 系统 的 土壤 保 
持 服务 、 产 水 服务 \ 碳 储存 服务 和 生境 质量 等 。 有 
关 生 态 系统 服务 的 研究 成 果 较 少 关注 矿区 这 类 特 
殊 的 生态 系统 , 且 现 有 的 少量 成 果 集 中 在 土地 类 型 
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转化 和 生态 脆弱 性 ”单个 时 间 节 点 的 定性 分 析 和 
静态 (或 截面 ) 研 究 ”, 缺 乏 长 时 间 序 列 矿区 生态 系 
统 服 务 价值 的 时 空格 局 分 析 ”。 矿 区 为 社会 经 济 
发 展 提供 了 大 量 资源 和 能 源 ,其 生态 系统 服务 研究 
正在 引起 学 术 界 的 关注 。 当 前 ,或 需 开 展 矿 区 生态 
系统 服务 价值 动态 变化 的 基础 研究 ,以便 更 好 地 服 
务 矿 区 资源 开发 利用 与 绿 水 青山 保护 的 协同 发 展 。 

平 朔 矿区 作为 中 国 最 大 的 露天 煤矿 开采 区 ,是 
内 煤炭 开采 区 生态 系统 变化 的 缩影 , 而且 地 处 黄 
土 高 原 边缘 和 北方 农 牧 交错 带 ,生态 环境 具有 脆弱 
性 和 过 渡 性 的 双重 特点 ,因此 在 平 朔 矿 区 开展 生态 
系统 服务 研究 具有 显著 的 代表 性 和 示范 性 。 本 研 
究 利 用 遥感 数据 和 土地 利用 数据 ,在 3S 技术 支撑 
下 ,采用 价值 当量 法 和 生态 贡献 度 等 方法 ,对 
1990 一 2020 年 平 朔 矿区 的 生态 系统 服务 进行 价值 
估算 ,定量 表征 1990 一 2020 年 平 朔 矿区 生态 系统 服 
务 价 值 的 时 空 变化 ,这 不 仅 可 以 为 黄土 高 原矿 区 开 
展 土地 复旦 和 生态 重建 提供 参考 依据 ,而 且 可 以 拓 
宽 生态 系统 服务 的 研究 范围 。 


1 数据 与 方法 


1.1 研究 区 概况 

本 研究 选取 黄土 高 原 典 型 矿区 一 一 平 朔 矿区 
为 研究 区 域 。 平 朔 矿 区 (39°18'~39°37’N,112°10’~ 
112s*28 中 ) 位 于 山西 省 北部 的 朔 州 市 境内 ,总 面积 约 
380 km。 气 候 类 型 为 中 温带 半 干 旱 大 陆 性 季风 气 
候 , 年 均 温 5$.0~7.5 % ,年 均 降 雨量 410~450 mm ,年 均 
蒸发 量 1750-2550 mm ,年 均 风 速 2.5~5.0 m s 99, 
地 貌 类 型 为 黄土 丘陵 ,黄土 分 布 广泛 ,再 加 上 主要 
土壤 类 型 为 栗 钙 土 和 栗 褐 土 ,土壤 有 机 质 含量 低 、 
结构 差 . 抗 蚀 能 力 弱 ,水 蚀 风 蚀 剧 烈 ,是 黄土 高 原 的 
典型 生态 脆弱 区 的 。 平 朔 矿 区 作为 中 国 五 大 露天 
煤矿 之 一 ,是 国内 煤炭 开采 时 间 最 长 .空间 跨度 最 
大 的 露天 矿区 。 矿 区 探 明 地 质 储 量 约 为 126x108 t, 
是 我 国 首 批 煤炭 国家 规划 的 19 个 矿区 之 一 ,区 内 有 
3 座 露 天 矿 ( 安 太保 露天 矿 、 东 露天 矿 、 安 家 岭 露 天 
矿 ) 和 3 座 千 万 吨 级 现代 化 井 工 矿 ( 井 工 一 矿 、 井 工 
二 矿 和 井 工 三 矿 )( 图 1)。 
1.2 数据 来 源 

本 研究 所 需 的 基础 数据 主要 包括 土地 利用 数据 、 
土壤 数据 社会 经 济 统计 数据 等 多 源 数据 集 ( 表 1)。 


图 1 研究 区 示意 
Fig. 1 Schematic diagram of the study area 


将 以 上 数据 利用 ArcGIS 10.2 工具 箱 中 的 样 条 
函数 捅 值 工具 (Spline) 和 数据 转换 工具 (Conversion 
Tools) 进 行 数 据 格式 上 的 均一 化 处 理 , 最 终 转换 为 
分 辨 率 30m 的 栅 格 数据 。 

1.3 研究 方法 

131 生态 系统 服务 价值 估算 方法 借鉴 Daily 和 
Costanza 等 评估 生态 系统 服务 价值 的 原理 和 方法 “3 , 
谢 高 地 团队 对 国内 200 位 生态 学 专家 进行 问卷 调 
A ,修订 了 耕地 估计 过 低 .湿地 估计 偏 高 等 问题 , 先 
后 于 2003 年 和 2008 年 制定 了 针对 中 国 陆地 的 生态 
系统 单位 面积 生态 服务 价值 当量 因子 表 天 2 。2015 
年 谢 高 地 团队 结合 文献 调研 .专家 知识 和 生物 量 等 
方法 ,修改 和 更 新 了 价值 当量 因子 表 ”, 这 是 国内 目 
前 最 科学 、 最 系统 的 价值 当量 因子 表 ”。 一 个 当量 因 
子 等 于 1 he 耕地 1a 所 产 粮食 的 市 场 价值 1/7”。 考 
虑 到 生态 系统 服务 具有 时 空 异 质 性 ,为 使 研究 结果 
更 切合 平 朔 矿 区 实际 ,本 研究 粮食 单产 量 采用 2020 
年 平 朔 矿 区 平均 粮食 单产 量 (5425.50 kg* hm ?) , 粮 
食 价 格 采 用 2020 年 全 国 玉米 ( 平 朔 矿 区 主要 粮食 
作物 是 玉米 ) 最 低 收购 价格 (1.20 元 "kg )。 经 计 
算 , 平 朔 矿区 生态 系统 服务 价值 的 一 个 当量 因子 为 
930.09 7G -hm^., RIE 4 E A TRA MIPEHS CIT 
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表 1 研究 数据 信息 


Tab.1 Information of study data 


数据 类 型 数据 要 求 及 处 理 数据 来 源 
土地 利用 数据 ”选取 1990 年 和 2000 年 Landsat TM 影像 .2010 年 和 2020 年 Landsat OLI 影 像 ,空间 分 辨 International Scientific Data Ser- 
率 均 为 30 m。 利 用 Erdas 9.1 软 件 的 多 波段 遥感 数据 融合 方法 ,对 各 期 影像 进行 (TM Vice Platform 
4-3-2 波 段 合成 ,OLL5-4-3 波 段 合成 ) 标 准 假 彩色 合成 ,二 次 多 项 式 校正 ,最 邻近 法 网址 为 :httpVydatamiroresdb en/ 
重 采 样 ,校正 误差 控制 在 1 个 像 元 以 内 中。 利用 ArcMap 9.3 软件 进行 土地 利用 的 人 机 
交互 判读 与 提取 , 解 译 显示 精度 控制 为 1:100000, 拓 扑 和 改 错 处 理 在 ArcInfo 9.3 环境 
下 完成 ,形成 图 形 数据 库 和 属性 数据 库 ”%。 为 了 提高 解 译 结果 的 准确 性 ,笔者 于 2020 
年 9 月 20 一 25 日 进行 野外 实地 调查 ,随机 抽取 不 同 地 类 各 20 个 样 点 ,进行 GPS 定 位 、 
相机 拍摄 景观 照片 ,对 照相 应 地 点 的 遥感 影像 ,现场 确定 地 类 ,与 解 译 结果 进行 对 比 
分 析 ,4 期 解 译 精度 都 在 92% 以 上 ,满足 研究 要 求 。 本 研究 土地 利用 划分 为 耕地 、 林 
地 草地 水域 .建设 用 地 及 未 利用 地 6 大 类 
土壤 数据 基于 世界 土壤 数据 库 (HWSD ) 的 中 国土 壤 数 据 集 (v1.1) ,分 辩 率 为 1 km Harmonized World Soil Database 
网 址 为 :http://www.fao.org/ 
社会 经 济 粮食 作物 种 类 面积 .产量 .价格 等 《 朔 州 市 统计 年 鉴 兴 全 国 农产品 
统计 数据 成 本 收益 资料 汇编 》 及 实地 调研 


当量 因子 值 与 自然 条 件 ,建设 用 地 借鉴 龙 精 华 等 
研究 成 果 , 可 以 得 到 平 朔 矿区 各 类 土地 生态 系统 服 
务 价 值 系数 ( 表 2) ,进而 计算 出 平 朔 矿区 土地 生态 
系统 服务 价值 。 计 算 公 式 为 

ESV= X (AXVC.) (1) 

ESV,- X (AxVC,) (2) 
式 中 :ESV ESV, 分 别 为 平 朔 矿区 土地 生态 系统 的 服 
务 价值 .第 y 项 服务 价值 (元 :a7);4, 为 平 朔 矿区 第 x 
种 地 类 面积 (hm) ;VC,、VC, 分 别 为 平 朔 矿区 第 x 种 
地 类 生态 系统 的 服务 价值 系数 .第 y 项 服务 价值 系 
数 ( 元 .hm?.a-)。 
1.3.0 生态 贡献 度 生态 贡献 度 表 示 在 某 一 时 间 段 
内 ,各 地 类 生态 系统 服务 价值 的 变化 量 对 平 朔 矿区 
总 服务 价值 变化 量 的 影响 ,可 以 揭示 影响 平 朔 矿区 
生态 系统 服务 价值 变化 的 主要 贡献 因子 和 敏感 因 


To HAARA”, 
CR, = 


u 


AESV, 
YjAESV 
il 


x 10096 (3) 


«| 

式 中 :CR, 为 第 ;类 生态 系统 在 :时 间 段 内 生态 贡献 
度 ;AESV; 为 第 i 类 生态 系统 在 1 时 间 段 内 的 生态 系 
统 服务 价值 变化 量 (元 "a ') ;n 为 地 类 (n=6)。 

1.3.3 敏感 性 分 析 引入 经 济 学 中 基于 弹性 系数 的 
敏感 性 指数 (CS) 定 量 描述 平 朔 矿区 生态 系统 服务 
价值 变化 对 价值 系数 的 依赖 程度 ,以 减少 研究 结 
的 不 确定 性 。 将 各 地 类 的 价值 系数 分 别 调整 50%， 
PEET KAH CS. RARA, 
(ESV' - ESV)/ESV 

(VC, - VCNC, 
式 中 :CS 为 敏感 性 指数 ;ESV ESV' 分 别 为 调整 前 、 


CS= (4) 


表 2 平 朔 矿区 各 类 土地 生态 系统 服务 价值 系数 


Tab.2 ESV coefficients of each land use type in Pingshuo mining area [J6*hm a 
一 级 类 型 二 级 类 型 耕地 林地 草地 水 域 建设 用 地 未 利用 地 
供给 服务 食物 生产 790.58 176.71 353.43 744.07 0.00 0.00 
原料 生产 372.04 399.94 520.85 213.92 0.00 0.00 
调节 服务 气体 调节 623.16 1131.43 1832.28 716.17 -831.90 167.42 
气候 调节 334.83 3934.28 4845.77 2129.91 0.00 502.25 
水 源 涵养 251.12 3115.81 3552.94 95092.40 -8550.70 27.90 
净化 环境 93.01 1190.52 1599.76 5161.99 -2087.80 93.01 
支持 服务 保持 土壤 957.99 1599.75 2232.22 864.98 4455.90 18.60 
生物 多 样 性 120.91 1460.24 2027.59 2371.73 0.00 18.60 
文化 服务 美学 景观 55.81 641.76 892.89 1757.87 0.00 93.01 
总 计 3599.45 13650.44 17857.73 109053.04 -7004.20 920.79 
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后 的 生态 系统 服务 价值 (元 :ai);VC VC 分 别 为 调整 
前 、 后 第 ;类 生态 系统 服务 价值 系数 (元 hm?.a-)。 
当 CS<1 时 ,表明 ESV 对 于 价值 系数 缺乏 弹性 ,计算 
结果 可 信 ; 当 CS>1 时 ,表明 ESV 对 于 价值 系数 富有 
弹性 ,计算 结果 准确 性 较 差 。CS 的 绝对 值 越 接近 0， 
准确 性 越 高 。 


2 结果 与 分 析 


2.1 平 朔 矿 区 土地 利用 的 时 空 变化 

通过 解 译 遥 感 影 像 ,得 到 平 朔 矿区 1990、2000、 
2010 年 和 2020 年 土地 利用 现状 图 (图 2) 和 土地 利 
用 统计 表 ( 表 3)。 

从 图 2 可 以 看 出 ,过 去 30 a 间 6 种 地 类 的 空间 分 
布 格局 ,耕地 始终 是 主要 地 类 ,上 且 空间 分 布 变化 不 
大 。 受 煤炭 露天 开采 的 扰动 ,2000 一 2010 年 安 太 保 
RAT 、 东 露天 矿 和 安家 岭 露天 矿 的 部 分 耕地 转化 
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为 建设 用 地 和 未 利用 地 ,从 中 部 出 现 并 逐渐 向 东 
部 、 东 北部 扩展 。 林 地 、 草 地 主要 分 布 在 井 工 一 矿 
和 井 工 三 矿 ,但 井 工 一 矿 的 林地 和 草地 转化 为 耕 
地 , 井 工 三 矿 的 林地 和 草地 逐渐 转化 为 建设 用 地 和 
未 利用 地 。 这 种 变化 格局 充分 体现 了 平 朔 矿区 从 
农业 为 主 向 农工 矿业 并 重 的 转变 ,使 得 本 区 成 为 具 
有 农工 矿 特 色 的 半 自 然 半 人 工 生 态 系统 结构 。 
2010—2020 年 平 朔 矿区 加 大 土地 复明 和 生态 修复 
力度 ,复明 区 逐渐 恢复 为 耕地 和 林 草 地 ,植被 覆盖 
率 达 到 90% 以 上 ,建设 用 地 面积 缩减 ,实现 了 从 匾 
山 废 坑 向 绿 水 青山 的 跨越 。 

从 变化 规模 来 看 ( 表 3),1990 一 2020 年 平 朔 矿 
区 的 耕地 面积 最 大 , 占 矿 区 总 面积 的 70% 以 上 。 随 
着 煤炭 开采 力度 的 加 大 和 城镇 建设 步伐 的 加 快 ， 
1990—2010 年 耕地 .草地 和 林地 面积 分 别 减少 了 
4081.29 hm^,1499.50 hm^,568.95 hm , 相 比 之 下 , 建 
设 用 地 .未 利用 地 和 水 域 均 处 于 增长 状态 ,最 为 显 
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Fig. 2 Present situation of land use in Pingshuo mining area from 1990 to 2020 


表 3 1990 一 2020 年 平 朔 矿区 土地 利用 变化 
Tab.3 Changes of land use in Pingshuo mining area from 1990 to 2020 


sn 1990 4E 2000 年 2010 年 2020 年 

面积 /hm 比例 /% 面积 /hm 比例 /% 面积 /hm 比例 /9p 面积 hm? 比例 /% 
耕地 31734.81 83.51 31234.12 82.19 27653.52 72.71 30612.09 80.56 
林地 1042.15 2.74 950.15 2.50 473.20 1.25 735.03 1.93 
草地 4103.62 10.80 3864.24 10.18 2604.12 6.85 3138.32 8.26 
水 域 2.12 0.01 5.13 0.01 63.25 0.17 15.04 0.04 
建设 用 地 930.06 2.45 1632.12 4.29 6735.56 17.73 2974.29 7.83 
未 利用 地 187.24 0.49 314.24 0.83 470.35 1.23 525.23 1.38 
合计 38000.00 100.00 38000.00 100.00 38000.00 100.00 38000.00 100.00 


著 的 是 建设 用 地 面积 由 1990 4F. 930.06 hm 增加 到 
2010 年 6735.56 hm , 净 增 加 5805.50 hm’, 2010 一 
2020 年 耕地 .草地 和 林地 面积 增加 ,而 建设 用 地 面 
职 减 少 (-3761.27 hm) 。 从 变化 率 来 看 ,过 去 30 a 
间 平 朔 矿区 建设 用 地 变化 幅度 最 大 , 占 比 由 2.45% 
上 升 到 7.83% ,增加 了 5.38% ; 耕地 和 草地 变化 幅度 
次 之 ,变化 率 分 别 为 -2.95% 、-2.54%; 林 地 、 未 利用 
地 和 水 域 变化 幅度 不 明显 。 

2.2 平 朔 矿区 生态 系统 服务 价值 的 时 间 变 化 

2.2.1 总 价值 的 时 间 变 化 随 着 煤炭 资源 大 规模 开 
发 ,1990 一 2020 年 平 朔 矿区 土地 生态 系统 服务 价值 
经 历 了 先 减少 后 增加 的 变化 过 程 , 由 1990 年 (资源 
开采 初期 ) 的 19562.43x10: 元 减少 到 2010 年 的 
11265.40x10' 7C ,再 增加 到 2020 4E. ÉY 15755.47x 10* 
元 ,变化 量 为 -3806.96x10' 元 ,变化 率 为 -19.46%( 表 
4). JEN .林地 .草地 和 建设 用 地 的 生态 系统 服务 
价值 与 平 朔 矿区 土地 生态 系统 服务 总 价值 变化 过 
程 相 似 ,变化 量 分 别 为 -404.12x10: 元 -419.23x10; 
J ,-1723.80x10* JÙ , -1431.82x 10* Jù ,变化 率 分 别 
为 -3.53% ,-29.4796 ,-23.5296 ,- 219.7996 ;水 域 和 未 
利用 地 的 生态 系统 服务 价值 变化 量 和 变化 率 均 为 
正 数 ,但 由 于 其 面积 和 变化 量 较 小 ,对 平 朔 矿区 土 
地 生态 系统 服务 总 价值 变化 影响 不 大 。 总 体 来 看 ， 
草地 面积 的 减少 和 建设 用 地 面积 的 增加 ,加 之 二 者 
生态 系统 服务 价值 系数 较 高 ,草地 和 建设 用 地 的 生 
态 系统 服务 价值 变化 量 较 大 (合计 -3155.62x107 
元 ) ,草地 和 建设 用 地 的 面积 变化 是 平 朔 矿区 生态 
系统 服务 总 价值 减少 的 主要 原因 。 
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从 分 时 段 来 看 ,1990 一 2000 年 平 朔 矿区 生态 系 
统 服务 总 价值 变化 量 为 -1180.50x10: 元 ,其 中 建设 用 
地 、 草 地 变化 量 较 大 ,分 别 为 -491.74x10’、-427.47x 
10! 元 ,变化 率 分 别 为 -75.49%、-5.83%。2000 一 2010 
年 总 价值 变化 量 为 -7116.53x10! 元 ,其 中 建设 用 地 、 
草地 变化 量 较 大 ,分 别 为 -3574.55x10’、-2250.29x10 
元 ,变化 率 分 别 为 -31.26% 、-32.61%。2010 一 2020 
年 总 价值 变化 量 为 4490.07x10' 元 ,其 中 建设 用 地 、 
耕地 变化 量 较 大 ,分 别 为 2634.47x10*、1064.92x10” 
元 ,变化 率 分 别 为 55.84% 、10.69%( 表 4)。 可 以 看 
出 ,在 3 个 研究 时 段 内 建设 用 地 和 草地 的 变化 量 最 
大 ,再 次 说 明 建 设 用 地 .草地 的 面积 变化 和 价值 量 
变化 对 平 朔 矿区 生态 系统 服务 价值 变化 造成 了 较 

影响 。 其 原因 是 , 煤 痰 资源 开采 改变 了 原 地 貌 ， 
草地 逐渐 被 工矿 建设 用 地 代替 ,导致 草地 退化 ,最 
终 造成 生物 多 样 性 减少 ,生态 系统 结构 改变 、 功 能 
退化 和 服务 价值 下 降 。 
2.2.2 单项 价值 的 时 间 变 化 1990 一 2020 年 平 朔 矿 
区 单项 生态 系统 服务 价值 变化 中 , 除 保持 土壤 价值 
增加 (11.88%) 外 ,其 他 各 项 生态 服务 价值 均 有 不 同 
程度 的 减少 ,其 中 水 源 涵养 .净化 环境 .生物 多 样 性 
价值 减少 幅度 较 大 ,分 别 减少 了 115.84% .69.929% 、 
18.29%( 表 5)。 分 时 段 变 化 中 ,1990 一 2000 年 平 朔 
矿区 单项 生态 系统 服务 价值 减少 最 多 的 是 水 源 涵 
FE .净化 环境 和 气候 调节 ,分别 减少 了 697.63x10* 
JÈ ,197.74x10* JÈ ,160.25x10*JÙ ; 2000—2010 年 单项 
价值 减少 最 多 的 是 水 源 涵养 净化 环境 和 和 气体 调 
节 ,分 别 减 少 了 4496.92x10: 元 1325.72x 10° JÈ 


表 4 1990 一 2020 年 平 朔 矿区 生态 系统 服务 价值 变化 
Tab.4 Changes of ESV in Pingshuo mining area from 1990 to 2020 


类 型 年 份 耕地 林地 草地 水 域 建设 用 地 未 利用 地 合计 

ESV/10' 元 1990 11422.79 1422.58 7328.13 23.12 -651.43 17.24 19562.43 
2000 11242.57 1296.99 6900.66 55.94 -1143.17 28.93 18381.93 

2010 9953.75 645.94 4650.37 689.76 -4717.72 43.31 11265.40 

2020 11018.67 1003.35 5604.33 164.02 -2083.25 48.36 15755.47 

ESV 变化 量 /10' 元 1990—2000 -180.72 -125.59 -427.47 32.82 -491.74 11.69 -1180.50 
2000—2010 -1288.82 -651.05 -2250.29 633.82 -3574.55 14.38 -7116.53 

2010—2020 1064.92 357.41 953.96 -525.74 2634.47 5.05 4490.07 

1990—2020 -404.12 -419.23 -1723.80 140.90 -1431.82 31.12 -3806.96 

ESV 变化 率 /% 1990—2000 -1.56 -8.83 -5.83 141.95 -75.49 67.80 —6.03 
2000—2010 -11.46 -50.19 -32.61 113.30 -31.26 49.71 -38.72 

2010—2020 10.69 55.33 20.51 -76.22 55.84 11.66 39.85 

1990—2020 -3.53 -29.47 -23.52 609.42 -219.79 180.51 -19.46 
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表 5 1990 一 2020 年 平 朔 矿区 单项 生态 系统 服务 价值 变化 
Tab. 5 Changes of individual ESV in Pingshuo mining area from 1990 to 2020 
供给 服务 调节 服务 支持 服务 文化 服务 
类 型 年 份 食物 ”原料 气体 “气候 ”水 源 ”净化 保持 ÆW 美学 合计 
生产 生产 调节 调节 涵养 环境 土壤 ”多样 性 景观 
ESV/10' 元 1990 2672.49 14362 . 277331 3470.96 1805.05 884.37 4538.83 1368.78 612.52 19562.43 
2000 2623.05 1401.41 2631.78 3310.7] 1107.42 686.63 4735.07 1301.71 584.15 18381.93 
2010 2291.34 1184.74 1706.02 2418.02 -3389.50 -639.09 6313.82 94734 432.1 11265.40 
2020 2545.16 — 1332.07 2331.31 2864.51 -285.96 265.96 507833 1118.33 505.76 15755.47 
ESV 变 化 1990—2000 -49.44 -34.71 -141.53  -160.25 -697.63  -197.74 196.24 — -67.07 -28.37 -1180.50 
量 /10 元 ^ 5999 2010 -331.71 -216.67 -925.76 -892.69 -4496.92 -1325.72 1578.75 -354.37 -151.44 -7116.53 
2010—2020 253.82 147.33 625.29 446.49 3103.54 905.05 -1235.49 170.99 73.05 4490.07 
1990—2020 -127.33 -104.05 -442.00 - 606.45 -2091.01 -618.41 539.50 -250.45 -106.76 -3806.96 
ESV 变 化 1990—2000 -1.85 -240 -5.10 -4.60 | -40.00 | -2240 4.30 . -490 -4.60 -6.04 
率 /% 2000—2010 -12.64 -15.46 -35.17 -26.96  -406.07 -193.07 33.34 -2722 -25.92 -38.71 
2010—2020 11.07 12.44 36.65 18.46 -91.56 -141.61 -19.56 18.04 16.88 39.85 
1990—2020 -4.76 -7.24 -15.93  -1747 -115.84 -69.92 11.88 -18.29 -17.43 -19.46 


925.76x10: 元 ;2010 一 2020 年 除了 土壤 保持 价值 外 ， 
其 他 单项 价值 均 有 所 增加 ,增加 最 多 的 是 水 源 涵 
JE .净化 环境 和 气体 调节 ,分 别 增加 了 3103.54x10* 
元 .905.05x10! 元 .625.29x10!: 元 。 总 体 来 看 ,调节 服 
务 价值 减少 驱使 平 朔 矿区 生态 系统 服务 总 价值 降 
低 , 而 保持 土壤 价值 的 变化 是 由 于 建设 用 地 减少 了 
水 土 流失 。 

从 平 朔 矿区 单项 生态 系统 服务 价值 的 构成 来 
看 ,过 去 30 a 间 各 单项 生态 服务 价值 居于 前 4 位 的 
是 :保持 土壤 气候 调节 食物 生产 和 气体 调节 ,四 
者 价值 之 和 占 各 年 期 土地 生态 系统 服务 总 价值 的 
70% 以 上 ,共同 维护 了 平 朔 矿区 生态 系统 服务 的 主 


要 功能 。 
2.3 平 朔 矿 区 生态 系统 服务 价值 的 空间 变化 

TE ArcGIS 10.2 软件 中 ,基于 机 格 统 计 了 1990, 
2000 , 2010 年 和 2020 年 平 朔 矿区 生态 系统 服务 价 
值 ,根据 服务 价值 高 低 分 为 5 类 ,反映 平 朔 矿 区 生态 
系统 服务 价值 的 空间 变化 (图 3)。 总 体 来 看 , 平 朔 
矿区 生态 系统 服务 价值 呈现 明显 的 南北 高 .中 部 低 
的 特征 ,变化 趋势 与 林地 、 草 地 和 建设 用 地 的 分 布 
基本 一 致 。 生 态 系统 价值 高 值 的 区 域 分 布 在 南部 
和 北部 部 分 地 区 ,最 高 值 位 于 井 工 一 矿 。 生 态 系统 
服务 价值 的 低 值 区 域 与 建设 用 地 、 未 利用 地 分 布 较 
为 一 致 ,主要 分 布 于 安家 岭 露 天 矿 中 部 、. 东 露天 矿 


(a) 1990 年 (b) 2000 年 (c) 2010 年 (d) 2020 年 N 
R= 
- SERE 

S 
图 例 
。 矿区 点 位 
COPT KHR 
生态 系统 服务 价值 /104 元 
7 -2014.51 


À E 2014.50-.18.59 
KERER 
ZeD EEE 18.60-748.20 


IH 748.21-1607.20 
3km wm 1607.21-9814.77 


图 3 1990 一 2020 年 平 朔 矿区 生态 系统 服务 价值 空间 变化 
Fig.3 Spatial changes of ESV in Pingshuo mining area from 1990 to 2020 


北部 和 安家 岭 露 天 人 矿 东南 部 。 生 态 系 统 价值 中 值 
分 布 范围 较 广 , 主 要 分 布 在 平 朔 矿区 的 中 部 和 中 南 
部 ,包括 安家 岭 露天 矿 南部 、 东 露天 矿 南部 、 井 工 二 
D 安家 岭 圳 天 矿 北 部 。 

从 分 时 有 段 来 看 ,1990 一 2000 年 平 朔 矿区 生态 系 
统 服 务 价值 的 空间 分 布 变 化 不 大 ,主要 原因 是 处 于 
煤炭 资源 开采 初期 ,对 生态 系统 结构 和 功能 的 影响 
尚未 体现 出 来 ;2000 一 2010 年 随 着 安 太 保 露 天 矿 持 
续 开 采 和 东 露 天 矿 建设 ,大 面积 耕地 转化 为 建设 用 
地 ,生态 系统 价值 低 值 区 域 在 两 矿 出 现 并 扩展 ; 
2010 一 2020 年 生态 系统 价值 低 值 区 域 呈 缩减 趋势 ， 
原因 是 矿区 大 力 实施 了 生态 环境 治理 和 土地 复 垦 
等 工作 ,部 分 建设 用 地 恢复 为 耕地 和 林 草 地 。 可 以 
预见 ,未 来 一 段 时 间 , 随 着 矿区 生态 环境 治理 持续 
推进 ,生态 系统 服务 价值 的 下 降 趋 势 将 得 到 扭转 。 
2.4 平 朔 矿区 的 生态 贡献 度 

从 1990—2020 年 平 朔 矿区 不 同 地 类 生态 系统 
服务 价值 变化 的 贡献 度 ( 表 6) 可 以 看 出 ,草地 、 建 设 
用 地 和 林地 的 生态 贡献 度 最 大 ,分 别 为 -41.52% 、 
-34.49% 、-10.09%。 平 朔 矿区 生态 贡献 度 总 和 为 
-91.71% ,过 去 30 a 间 生态 系统 服务 价值 减少 了 
3806.96x10' 元 ,降幅 为 19.46%。 导 致 生态 系 统 服务 
价值 损失 的 地 类 有 建设 用 地 .草地 .林地 和 耕地 , 生 
态 贡 献 度 之 和 为 -95.84% ,建设 用 地 .草地 和 林地 是 
主要 的 负 向 贡献 因子 ;导致 生态 系统 服务 价值 增益 
的 地 类 有 水 域 和 未 利用 地 ,生态 贡献 度 之 和 为 
4.13% ,水域 是 主要 的 正 向 贡献 因子 。 建 设 用 地 、 草 
地 和 林地 是 平 朔 矿区 生态 系统 服务 价值 变化 的 主 
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要 贡献 因子 和 敏感 因子 。 分 析 3 个 研究 时 段 中 不 同 
地 类 的 生态 贡献 度 ,也 可 得 出 同样 的 结论 。 
2.5 平 朔 矿区 的 敏感 性 分 析 

从 敏感 性 指数 计算 结果 可 以 看 出 ( 表 7) ,价值 系 
数 调整 50% 后 , 平 朔 矿区 不 同 地 类 的 敏感 性 指数 区 
别 较 大 ,但 同一 地 类 不 同年 份 之 间 差 别 不 大 , 且 各 时 
期 敏感 性 指数 均 小 于 1。1990 一 2020 FET KA 
地 的 敏感 性 指数 最 高 ,其 范围 为 0.584~0.884, 即 当 
价值 系数 增加 10% 时 ,生态 系统 服务 总 价值 增加 
5.84%~8.84% ;未 利用 地 敏感 性 指数 最 低 ,接近 0; 当 
未 利用 地 的 价值 系数 增加 10% 时 ,生态 系统 服务 总 
价值 最 大 增加 0.04% ,波动 较 小 。 由 此 得 出 , 平 朔 矿 
区 生态 系统 服务 总 价值 对 价值 系数 缺乏 弹性 ,本 研 
究 所 采用 的 价值 系数 符合 平 朔 矿 区 的 实际 情况 , 研 
究 结 果 可 信 。 
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本 研究 基于 遥感 数据 和 土地 利用 数据 ,采用 3S 
技术 得 到 1990—2020 年 平 朔 矿区 土地 总 面积 不 变 
的 情况 下 ,各 地 类 面积 发 生 了 较 大 变化 (耕地 等 减 
少 和 建设 用 地 等 增加 )。 在 国家 政策 的 引导 下 , 平 
朔 矿 区 的 露天 煤矿 相继 投产 ( 安 太保 露天 矿 1987 年 
投产 ,安家 岭 露天 矿 2005 年 投产 , 东 露 天 矿 2016 年 
投产 ) ,原煤 产量 由 1990 年 的 1250x10t 增 加 到 2020 
ERY 8100x10* t, 人口 同 期 增加 了 7.20x10* 人 (https:// 
pingshuo.chinacoal.com/) 。 为 了 满足 煤炭 业 的 发 展 
和 居民 的 需求 ,当地 政府 增加 建设 用 地 用 来 建设 厂 


RO 1990 一 2020 年 平 朔 矿 区 不 同 地 类 生态 系统 服务 价值 变化 的 贡献 度 


Tab.6 Contribution rates of changes in ESV of different land use in Pingshuo mining area from 1990 to 2020 


/% 


年 份 耕地 林地 草地 水 域 建设 用 地 未 利用 地 
1990—2000 -14.23 —9.89 -33.66 2.58 -38.72 0.92 
2000—2010 -15.31 -1.73 -26.74 7:53 -42.48 0.17 
2010—2020 19.22 6.45 17.21 一 9.48 47.54 0.09 
1990 一 2020 -9.74 一 10.09 -41.52 3.39 -34.49 0.74 

表 7 1990—2020 年 平 朔 矿 区 生态 系统 服务 价值 敏感 性 指数 
Tab.7 Sensitive coefficients of ESV in Pingshuo mining area from 1990 to 2020 

年 份 耕地 林地 草地 水 域 建设 用 地 未 利用 地 

1990 0.584 0.073 0.375 0.001 0.033 0.001 

2000 0.612 0.071 0.375 0.003 0.062 0.002 

2010 0.884 0.057 0.413 0.061 0.419 0.004 

2020 0.699 0.064 0.356 0.010 0.132 0.003 
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矿 、 交 通道 路 和 生活 设施 等 。 这 种 土地 利用 变化 格 
局 与 朱 九 龙 ”\ 龙 精华 等 弛 对 2005 一 2012 年 河南 义 
马 煤 矿 矿 区 .1993 一 2013 年 黑龙 江 省 稚 岗 矿区 土地 
利用 变化 研究 成 果 一 致 。 

无 论 Costanza 等 5 提出 的 生态 系统 服务 价值 系 
数 法 ,还 是 谢 高 地 等 ”修改 和 更 新 后 的 价值 当量 因 
子 表 , 对 大 尺度 生态 系统 服务 研究 较为 可 信 , 但 在 小 
尺度 研究 中 由 于 各 种 自然 因素 和 人 文 要 素 差 异 较 
小 ,研究 结果 存疑 较 多 。 近 年 来 ,不 同学 者 ”* 采 
用 修正 后 的 单位 面积 价值 当量 法 对 不 同 区 域 进 行 
生态 系统 服务 价值 评估 且 得 到 了 较为 准确 的 评估 
结果 。 本 研究 以 价值 当量 法 为 基础 ,结合 平 朔 矿区 
实际 情况 ,利用 1990 一 2020 年 粮食 价值 对 当量 进行 
修正 ,与 生态 系统 服务 价值 的 静态 估算 方法 2” 相 
比 ,计算 结果 更 具有 准确 性 全面 性 和 分 析 性 ,得 到 
1990—2020 年 平 朔 矿区 生态 系统 服务 价值 呈现 先 
下 降 后 上 升 趋势 ,空间 上 呈现 南北 高 .中 部 低 的 分 
布 特征 ,这 种 矿区 生态 系统 服务 的 时 空格 局 具有 鲜 
明 的 地 域 性 和 生态 治理 的 必要 性 ”23 。 

矿区 生态 系统 包含 多 个 子 系统 .多 种 服务 , 且 
各 服务 之 间 相 互联 系 、 相 互 影响 ,任何 一 种 生态 系 
统 服务 的 变化 都 会 引起 其 他 服务 或 正 或 负 的 变 
化 。 因 此 ,基于 本 研究 结果 ,未 来 将 采用 更 科学 的 
数据 和 模型 深入 挖掘 平 朔 矿区 土地 利用 变化 和 生 
态 系统 服务 变化 的 驱动 机 制 ,为 矿区 复 垦 与 生态 重 
建 及 管理 政策 的 制定 提供 科学 依据 。 


4 结 论 

(1) 1990—2020 年 平 朔 矿区 的 耕地 .草地 和 林 
地 面积 呈现 缩减 趋势 ,分别 减少 了 1122.72 hm’, 
965.30 hm, 307.12 hm'; 建设 用 地 面积 扩展 最 大 
(2044.23 hm ) 。 

(2) 弥补 了 生态 系统 服务 研究 较 少 关注 矿区 
生态 系统 的 不 足 , 分 析 了 长 时 间 序 列 (30 a) 矿 区 
生态 系统 服务 价值 的 时 空格 局 。1990 一 2020 年 平 
朔 矿 区 生态 系统 服务 总 价值 呈 先 下 降 后 上 升 态 
势 , 由 1990 年 的 19562.43x10! 元 减少 到 2010 年 的 
11265.40x10' 7C ,再 增加 到 2020 4E. ÉY 15755.47x 10* 
元 。 平 朔 矿区 除 保 持 土壤 价值 增加 外 ,其 他 各 单项 
生态 服务 价值 均 有 不 同 程度 的 减少 。1990 一 2020 
年 平 朔 矿区 生态 系统 服务 价值 空间 分 布 呈现 南北 


高 .中 部 低 的 特征 ,空间 分 布 与 林地 .草地 和 建设 用 
地 的 分 布 基本 一 致 。 

(3) 1990—2020 年 平 朔 矿区 草地 建设 用 地 和 
林地 的 生态 贡献 度 最 大 ,分 别 为 -41.52% , 34.4996 
和 -10.09% ,三 者 是 主要 的 贡献 因子 和 敏感 因子 。 
1990—2020 年 平 朔 矿 区 的 敏感 性 指数 值 均 小 于 1, 
人 研究 结果 可 信 , 表 明基 于 当量 因子 修正 方法 计算 获 


得 的 生态 系统 服务 价值 系数 符合 平 朔 矿 区 的 实际 
情况 。 
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Spatiotemporal changes of the ecosystem service value for mining area in 
Loess Plateau: A case of Pingshuo mining area 


MENG Xianwen", CAO Jun, XUE Zhanjin 


(1. Taiyuan Tourism College, Taiyuan 030032, Shanxi, China; 2. Institute of Resources and Environmental Engineering, Shanxi 
University, Taiyuan 030006, Shanxi, China; 3. Institute of Ecological Civilization and Green Development, Guangdong Provincial 


Academy of Environmental Science, Guangzhou. 510000, Guangdong, China) 


Abstract: Spatiotemporal changes in ecosystem service value are crucial for protecting the ecological environ- 
ment and enhancing the ecosystem functionality within mining areas. Based on remote sensing and land-use data, 
ecosystem service value (ESV) was quantitatively analyzed in the Pingshuo mining area of the Loess Plateau of 
China from 1990 to 2020 using the 3S (GIS, RS, GPS) technique, a standard equivalent method, and ecological 
contribution rate. The results showed the following facts: (1) Shifts in land-use type from 1990 to 2020, notably 
the most substantial decrease occurring in farmland (1122.72 hm’) and a corresponding increase in built-up land 
(2044.23 hmô) to its highest. (2) A decline in total ESV from 19562.43x10" Yuan in 1990 to 11265.40x10* Yuan in 
2010, primarily attributed to reduced ESV in built- up land and grassland. However, a subsequent increase to 
15755.47x 10° Yuan in 2020 was observed, which was attributed to the implementation of ecological programs. 
(3) Varied decreases in individual ESV components, with water retention, purified environment, and biodiversity 
values decreased by 115.8495, 69.9296, and 18.29% , respectively, over 1990—2020, excluding soil retention, 
which remained constant. (4) ESV spatial distribution in the Pingshuo mining area, exhibiting higher values in the 
north/south and lower in the middle. (5) Ecological contribution rates of grassland, built-up land, and forestland 
in 1990—2020, standing at —41.5296, -34.49%, and —10.0995, respectively. These factors emerged as the main 
contributors and sensitive elements in the Pingshuo mining area. 


Key words: ecosystem service value; spatiotemporal change; Pingshuo mining area; Loess Plateau 


